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1 Introduction 



Over the last decade or two there has been an increasing interest in efficient quan- 
tization procedures for simple (having only few freedoms) but nonlinear (possibly 
chaotic) Hamiltonian systems. Here I consider bound and autonomous Hamiltonian 
systems with / freedoms. Directly solving the time-independent Schrodinger equa- 
tion in /—dim configuration space (CS) or the equivalent eigenvalue problem for the 
Hamiltonian matrix in an appropriate basis is the first but certainly not the best 
idea. The question appeared whether there exists a quantum analogue of surface 
of section (SOS) reduction of classical dynamics || which reduces smooth bound and 
autonomous Hamiltonian dynamics over 2/-dim phase space to a discrete Poincare 
mapping over only (2f — 2)-dim SOS. 

In the case of quantum billiard systems in two dimensions (/ = 2) we have 
the so called boundary integral method which reduces 2-dim Schrodinger equation 
to 1-dim integral equation. Its kernel can be interpreted as the quantum bounce 
map which is a special case of Poincare mapping. Smilansky and coworkers |2|, [14}| 
have developed more general scattering approach for quantization of billiards. They 
construct exact quantum Poincare mapping for 2-dim billiards with respect to the 
arbitrary line of section as the product of the two scattering matrices of the two 
opened billiards. These methods are typically much more efficient than the direct 
diagonalization, since the dimension of the matrices that they use is typically of the 
order of square root of the dimension of the original Hamiltonian matrix. On the 
other hand Bogomolny succeeded to construct approximate semiclassical Poincare 
mapping with respect to arbitrary configurational surface of section for arbitrary 
autonomous Hamiltonian. In this paper I present the generalization of the scattering 
approach for quantization of almost arbitrary bound Hamiltonians and show that it 
reduces to the Bogomolny 's theory in the semiclassical limit % — > 0. 

In section 2 I construct the quantum Poincare mapping and prove that the 
eigenenergies of the original Hamiltonian correspond to the fixed points of quantum 
Poincare mapping. I prove also theoretically perhaps even more interesting SOS 
decomposition of the resolvent of the Hamiltonian. Here I also study the semiclas- 
sical limit of newly defined propagators and explicitly calculate the leading order 
and next-to-leading order terms while I explain how higher order corrections (in 
powers of h) can be systematically obtained. The symmetries of quantum Poincare 
mapping are discussed and it is explained how the SOS quantization condition can 
be very efficiently used in practical calculations, especially for the generic class of 
the so called semi- separable systems. In section 3 I formulate an abstract quantum 
SOS method which can be applied to arbitrary boundary value differential equa- 
tion problems. Then I apply an abstract theory to the case of energy-dependent 
Schrodinger equation of section 2 and more general cases of nonrelativistic or even 
relativistic systems (described by the Dirac equation) coupled to arbitrary external 
gauge fields. In section 4 the method is generalized to the case of non-simply but 
multiply connected CSOS. In section 5 I discuss the meaning and applicability of 
the new results and conclude. 

Some preliminary results of this project have already been reported || [H|. Re- 
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cently I have been informed that one of the results of this paper, namely the SOS 
quantization condition for 2-dim Hamiltonian systems of the standard type, has also 
been obtained independently and subsequently by other authors [T^ . 



2 Surface of Section Quantization 
2.1 Notation 

The basic results of this paper are most beautifully and compactly written in terms 
of some new physical quantities whose mathematical definitions and notation are 
described in this subsection. 

We study autonomous and bound (at least in the energy region of our concern) 
Hamiltonian systems with few, say /, freedoms, living in an /—dim configuration 
space (CS) C. One should also provide a smooth (/ — l)-dim submanifold of CS C 
which shall be called configurational surface of section (CSOS) and denoted by 
Sq. In this section we consider only the case of simply- connected CSOS whereas in 
section 4 we study the case of more general multiply- connected CSOS. We choose 
the coordinates in CS, q = (x,y) G C in such a way that the CSOS is given by 
a simple constraint y = 0, or Sq = (5,0). These coordinates need not be global, 
i.e. they need not cover the whole CS, but they should cover the open set which 
includes the whole CSOS S . This means that every point in S should be uniquely 
represented by CSOS coordinates x 6 S which may be more general than Euclidean 
coordinates TV~ l (e.g. (/ — 1) dim sphere S^ 1 ). In this section we shall assume 
that So is an orientable manifold so that it cuts the CS C in two pieces which will be 
referred to as upper and lower and denoted by the value of the binary index a =|, j, 
(see figure 1). In arithmetic expressions the arrows will have the following values 
t= +1,|= —1. My approach presented in this section applies to quite general class 
of bound Hamiltonians whose kinetic energy is quadratic at least perpendicularly to 
CSOS 

H=^-pl + H'(p x ,x,y). (1) 

In next sections we generalize this class to include Hamiltonians having coordinate- 
dependent mass (which arise in the curvilinear coordinates which must be used in 
case of non-flat CSOS) and/or terms linear in p y (which appear e.g. due to the 
presence of magnetic field). 

In quantum mechanics, the observables are represented by self-adjoint operators 
in a Hilbert space Ti of complex- valued functions ^{q) over the CS C which obey 
boundary conditions ty{dC) = (] and have finite L 2 norm f c dq\^(q)\ 2 < oo. We 
shall use the Dirac's notation. Pure state of a physical system is represented by 
a vector — ket which can be expanded in a convenient complete set of basis 

2 The more general case of (2/ — 2)-dim SOS in 2/-dim phase space which is not perpendicular 
to CS cannot be treated within the present approach except in the cases where one can change the 
phase space coordinates by means of an appropriate canonical transformation. 

3 In case of infinite CS (e.g. Euclidean space C = W) we assume that points at infinity belong 
to the boundary dC. 
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vectors, e.g. position eigenvectors \q) = \x,y), = J e dq\q)(q\^) = J c dq^(q)\q) 
(in a symbolic sense, since \q) are not proper vectors, but such expansions are 
still meaningful iff ^{q) = (q\^) is square integrable i.e. L 2 (C)-function). Every 
ket l^) G Ti has a corresponding vector from the dual Hilbert space 7i', that is 
bra (^\ G 7i', (^\q) = (q\^)*. We shall use mathematical accent "to denote linear 
operators over the Hilbert space Ti. Operators of SOS coordinates x and p x , defined 

(x,y\x\ty) = x^(x,y), 
(x,y\Px\®) = -tfi>d x V(x,y), 

can be viewed also as acting on functions ip{x) of x only and therefore operat- 
ing in some other, much smaller Hilbert space of square-integrable complex- valued 
functions over a CSOS S 

{x\x\ip} = xtp(x), 
{^iPxl^} = -ihd x ip(x). 

Vectors in such SOS-Hilbert space, denoted by C, will be written as \ip} and linear 
operators over £ will wear mathematical accent " like restricted position x and mo- 
mentum p x . Eigenvectors \x} of SOS-position operator x provide a useful complete 
set of basis vectors of C. The quantum Hamiltonian can be written as 

H = -^-d 2 y + H'(y), H'(y) = H'(-ihd x ,x,y). (2) 

The eigenstates of the inside-CSOS Hamiltonian H'(0) = H'(0)\c restricted to the 
SOS-Hilbert space £, |n} G C 

F'(0)|n} = £>}, (3) 

which are called SOS-eigenmodes, provide a useful (countable n = 1,2,...) com- 
plete and orthogonal basis for C since H'(0) is a self-adjoint operator with discrete 
spectrum when its domain is restricted to C. 

The major problem of bound quantum dynamics is to determine the eigenener- 
gies E for which the Schrddinger equation 

(x,y\H\**(E)) = E* a (x,y,E) (4) 

has nontrivial normalizable solutions — eigenf unctions ty(x,y,E). 

2.2 Scattering formulation 

In this subsection I will introduce our basic tools using the powerful quantum me- 
chanical time-independent multi-channel scattering theory || [15]]. 

4 In case of non-euclidean SOS they should be replaced by the generators of the corresponding- 
Lie algebra. 
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To connect bound Hamiltonian dynamics and scattering theory one should make 
the following very important step. Cut one part of CS off along CSOS and attach 
a semi-infinite separable (flat along the y — axis) waveguide instead (see figure 1). 
Thus we introduce two scattering Hamiltonians 

~ f -(h 2 /2m)d 2 + H'(y); ay>0, 
° \ -(h 2 /2m)d 2 + H'(0); ay<0. 1 ' 

Every wavefunction inside the waveguide (ay < 0) at energy E can be separated as 
the superposition of products of a bound state (SOS-eigenmode n) in x direction 
and free motion in y— direction, 

{x\n}e ±iK{E)y , 

with the corresponding wavenumber determined by the energy difference E — E' n 
available for the motion perpendicular to the CSOS 



2m 

k n (E) = ] j w (E-E' n ). 

For any value of energy E, there is typically a finite number of the so called open 
or propagating SOS-eigenmodes — channels with real wavenumbers for which E' n < 
E, and an infinitely many closed channels with imaginary wavenumbers for which 
E' n > E. Scattering wavefunction ^ a (x, y, E) at a given energy E (or complex 
conjugated wavefunction ty* (x, y, E*) = (^ ff (x,y, E*))*) satisfying the Schrodinger 
equation H a \ty a (E)) = E\^/ cr (E)) can be uniquely parametrized by vector from 
the SOS-Hilbert space C (or by vector from the dual SOS-Hilbert space £). 
6 C essentially parametrize the incoming waves 



K(x,y,E) = ^{ 8 | n > / -^L— e ^(^{n|V>} 

Y in k n (E) 

which uniquely determine the whole scattering wavefunction. Therefore the wave 
operators can be defined, namely Q' a (E) which map from £ to 7i (or P^{E) which 
map from TC to C) and whose kernels are given by the scattering wavefunctions (or 
their complex conjugates) 

(q\Q' a (E)m = 9*(q,Ef), (6) 
mK{E)\q) = K(q,E*). (7) 

On the a— side of CS (ay > 0) the scattering wavefunction satisfies ordinary Schro- 
dinger equation (f|) whereas in the waveguide (ay < 0) it is a superposition of 
incoming and scattered waves 





-im 


h 




-im 


h 


V 


-im 



^a(x, y, E) = ^l^\n}k- 1/2 (E) [e^^Sra + e"^)^] {m = 

n,l 

'-{x\k- l ' 2 (E) \e iaR ^y + e~ laIi{E)y f a (E)} \tJ>}, (8) 



%(x,y,E*) = ^— ^{^*| \e iaR{E)y + f a (E)e- laR{E)y } R- 1/2 (E)\x}. (9) 
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For the sake of compact notation we have introduced the wavenumber operator 



2m 

K(E) = J2k n (E)\n}{n\ = - H'(0)). (10) 

T^(E) is the generalized scattering matrix since it includes also closed - - non- 
propagating modes and T a (E) is the corresponding scattering operator over £ 

T<r{E) =Y, T ni\n}{l\- (11) 

n,l 

Here I have to make three important notes: 

• Conjugated energy E* is used in argument of complex conjugated wavefunction 
(0) in order to make all the relevant operators, e.g. P' a (E), complex analytic 
functions of E rather than E*. 

• The SOS-states \ip*} and are generally different. 

• The equation @ is non-trivial and it does not follow from eq. (||) but it is a 

consequence of the Hermitian symmetry of the scattering Hamiltonians H a as 
will be shown in the next paragraph. 

Let us now consider the resolvents of the scattering Hamiltonians @ with out- 
going boundary conditions 

GJE) = — / dte i{E - H ° )t/h = (E-H a + tO)- 1 . (12) 
in Jo 

It is convenient to introduce a hybrid representation of these scattering Green func- 
tions denoted by G a (y,y',E) G £ (being a matrix element in y— variable and an 
operator in cc-variable) defined as 

{x\G a (y,y',E)\x'} = (x,y\G a {E)\x',y'). (13) 

Inside the waveguide (ay < 0, ay' < 0) these hybrid Green functions satisfy the 
following "free motion" Schrodinger equations in both arguments 

v 2 777 

d 2 G a (y,y',E) + K 2 (E)G a (y,y',E) = -^6(y-y'), (14) 

d 2 G a (y,y',E) + G a (y,y',E)K 2 (E) = -^5{y-y'). (15) 

The general solution of this linear system (in the waveguide) is given by the sum of 
particular — "free motion" solution 

GUv,y',E) = ^K^ 2 {E)e^\y-y'\k-^ 2 {E) 

in, 

and general solution of the homogeneous system satisfying outgoing boundary con- 
ditions 

G a (y, y', E) - G bee (y, y', E) = ^fr^We**™* 'Ae'™*™ 'k^E). (16) 

in 



Their sum {x\G a (y, 0, E) satisfies the Schrodinger equation ([|) on the a— side {ay > 
0), so comparing it locally, at ay = +0, with wavefunctions (H) yields 



V*{x,y,E) 



V- 



■im 



{x\G a {y^,E)k l l 2 {E)\i,}, ay>0, 



(17) 



and determines the free operator valued parameter, A = T a (E). Thus the waveguide 
expression for the hybrid scattering Green function reads (ay < 0, ay' < 0) 



G a (y,y',E) = ™K- X I\E) UWI^'I + e ^*(^vf ff (£) e -wW 
in L 



R-^iE). 



(18) 

Since G a (0,y, E)\x} satisfies conjugated Schrodinger equation, there exist SOS- 
states such that 



n 



im 



-.{r\K 1/2 (E)G a (0,y,E)\x}, ay > 0, 



(19) 



and equation <M) follows. 



2.3 SOS energy quantization 

Now I shall formulate an exact energy quantization condition for the original Hamil- 
tonian matrix H solely in terms of the scattering operators T a (E). 

Theorem la: Every energy E for which the operator 1 — T^(E)T^(E) (where the 
order of the arrows may be reversed) is singular is 

• either eigenenergy of the original Hamiltonian H , 

• or it is a threshold energy for opening of a new channel, 

• or both. 

More precisely: The dimensions of the left and right null-space of an operator 1 — 
Ti(E)T^(E) are the same 

d T (E) = dimker(l - f^f^E)) = dimker(l - f^f^E)^ (20) 

and the following inequality for dx{E) in terms of dimension of the null space of 
operator E — H, dn{E) = dimker(.E — H), and dimension of the null space of 
operator K 2 (E) } d K (E) = dimker K 2 (E), holds 

max{d H (E), d K (E)} < d T (E) < d H (E) + d K (E). (21) 



Proof: Let driEo) SOS-states | f n} G C, n = 1, . . . , dr{Eo) span the null space of 
l-T^W) 

(l-f [ {E Q )f^E ))\]n} = 0. (22) 
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Then one may define another set of (It(Eq) SOS-states | j n} G C by the prescription 

| |n} = f T (^)| tn}, (23) 
in terms of which the equation (|22| ) may be rewritten as a relation symmetric to 

(H 

|Tn}=f i (^)Un}. (24) 

Each of these vectors | f n} lies either in the null space or in the image of K 2 (Eq), 
since 

£ = ker K 2 (E )®K 2 (E )C, 

where K 2 (E)C = kerK 2 {E)- L since K 2 (E) is self-adjoint. Let the first rrix vectors 
| m}, m = l,...,rriK lie in ker K 2 (E). In order to make sure that scattering 
wavefunctions (§J^) have regular limit E — > i?o, since K~ x l 2 {Ei) is becoming singular 
if ker^T 2 (£'o) 7^ 0> one should demand 

(1 + T„(£ )M = 0, {0|(1 + f a (E )) = 0, (25) 

for any \<f)} G keri^ 2 (i? ), so ker i^ 2 (i? ) is invariant under T a (E ) and T^(E ) 

f a (E ) ker K 2 (E ) = f^(E ) ker K 2 (E ) = ker K 2 (E ). 

From (^5|) we also see that keri^ 2 (i?o) C ker(l — T^(E )T^(E )), so | | m } span 
the entire space keri^ 2 (£' ), and so cIk(E ) = m K < cIt(E ). Therefore the image 
K 2 (E)C is also invariant under T a (E ) and T^(E ), so the counterparts | j 1} of the 
remaining dr(Eo) — cIk(Eq) SOS-states | | I}, I — cIk(E ) + 1, . . . , cIt(E ) from the 
image K 2 [E)L also lie in the image K 2 (E)C In the image one can define the inverse 
of K 2 (E ) and the inverse of its fourth root, namely K~ 1 I 2 (Eq). Using (^,^|) one 
can write 

(l+TV(JSb)) | TO = (l + f i {E Q ))\{l} 1 

(l-T T (£ ))|T0 = - (l - to)) I I «}■ 

which can be rewritten using the values of the wavefunctions and their normal 
derivatives on the CSOS (from @) 

(^OIQ'^M = ^^{^^-^(^(l + f^))^}, (26) 
d y {x,y\Q'Am}\y=o = a^pixlK^tl-UE))^} (27) 
as the continuity of the wavefunctions and their normal derivatives on the CSOS 

(x,0\Q\(E )\ TO = (x, 0\Q{(E )\ 11} 
d y (x,y\Q'JE )\U}\y=o = d y (x,y\QUE )\ll}\ y=0 . 
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which are built up from the pairs of scattering wavefunctions 



(Ih(Eo), the maximal number of such linearly independent eigenfunctions ^i(x,y) 
is at least d T (E ) — d K (E ) since due to completeness of SOS-eigenstates |n} the 
mapping Q' a (E ) is injective on the image K 2 (E )C But this number, cIh(Eo), can 
be larger than dr(Eo) — (Ik{Eq) (but not larger than dx(E )) since there may be 
some states from the null space of K 2 (E ) for which the limits (^,^) of the upper 
(a =f) and lower (a =|) part accidentally match. Therefore we have proved an 
inequality 

Analogously, one can show that the d' T (E Q ) basis vectors | j n*} G C, n = 
1, . . . , d' T (Eo) of the left null space 

{in^l-f^)^)) =0. (29) 
are mapped onto conjugated basis of eigenfunctions under the propagator P^(E ) 

*» (a! ' y) -\ Un1^)|«,y>; y < 0, (3 ° } 

where the counterparts | | n*} are again defined as 

{Tn*| = Un*|f ; (^). (31) 

Generally, these conjugated wavefunctions are continuous and differentiate on the 
CSOS and are thus eigenfunctions of the Hamiltonian H if \an*} G K 2 (E )C, oth- 
erwise if \an*} G ker K 2 (E ) the continouity can be accidentally satisfied in the 
limit E — > Eq if both contributions (upper and lower) coincide. This happens 
when the corresponding limits for |cm} (^,[27]) coincide since the two cases differ 
only by a complex conjugation. So, the dimensions of left and right null space of 
1 — T^(Eq)T^(E ) should be the same dT(E ) = d' T (E ). 



The operator T'(E) = T^(E)T-^(E) will be called quantum Poincare mapping and it 
is the product of the two quantized Poincare scattering mappings. We have proved 
an extremely efficient quantization condition (as we shall show later), namely: The 
energies where the quantum Poincare mapping has fixed points (eigenvalue 1) are 
either: (i) eigenenergies of the Hamiltonian H or (ii) thresholds for opening of new 
channels E' n (which are already known as a solution of (|3|) as a prerequisite of the 
method). 



2.4 SOS decomposition of the resolvent of the Hamiltonian 

The kernels of the scattering propagators {x\T a (E)\x'} will henceforth be called 
CSOS-CSOS propagators. Then we define also: (i) The linear operator Q a (E) from 
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£ to 7i and the linear operator P<j{E) from 7i to £ with the kernels 

(x,y\Q«(E)m = \ frvfoWM OVZO (32) 

^ u, try <c u, 

{^(£)|a^> = {{fl^l^v); *v^o, (33) 

which are called CSOS-CS and CS-CSOS propagators respectively, and (ii) a linear 
operator Gq(E) over 7i with the kernel 

f (aj^lGft^laj'^'); 2/ >0, 2/>0, 
(*,y|6 (£) = ^ (£)!*', y'); 2/ < 0, y' < 0, (34) 

I 0; yy' < 0. 

which is called CS-CS propagator (without crossing the CSOS in between). 

Theorem 2a: The energy- dependent quantum propagator (i.e. the resolvent of the 
Hamiltonian) G(E) = (E — H)^ 1 can be decomposed in terms of the CS-CS prop- 
agator — with no intersection with the CSOS So — Gq{E), CS-CSOS propagator 
P a (E), CSOS-CS propagator Q a (E) , and CSOS-CSOS propagator T a (E) 

G(E) = G (E) + £Q CT (£)(1 -f^{E)f a {E))- l P^{E) 

(J 

+ ^^(^(l-T^^f^E))- 1 ^^)^^). (35) 



Quantities (q\G (E)\q') , (q\Q a (E)\x'}, {x\P a (E)\q') , {x\f a (E)\x'} should be in- 
terpreted as the probability amplitudes to propagate through the cx-side of CS from 
point q' in CS / x' on CSOS to point q in CS / x on CSOS at energy E and with- 
out crossing CSOS in between. Then this decomposition formula can be intuitively 
understood by expanding the operator (1 — T_ cr (_E , )T cr (£'))~ 1 in a geometric series 
and then using the basic postulates of quantum mechanics about summation of the 
probability amplitudes of alternative events (different number of crossings of CSOS) 
and multiplication of the probability amplitudes of consecutive events (sequential 
crossings of CSOS) [|J, since the system which propagates from point q i to point q^ 
in CS along continuous path can cross the CSOS arbitrarily many times. (In fact, 
the number of crossings is even if q { and q^ lie on the same side of CSOS and odd 
otherwise.) Two versions of the proof of this formula are given in || |TD[ while in 
this paper the proof will be given for more general cases which include the present 
one in the following two sections. 

2.5 Semiclassical limit 

In order to find explicit leading order semiclassical expressions for the CSOS/CS- 
CSOS/CS propagators it is convenient to express them first in terms of the scattering 
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Green functions in the hybrid representation ( jl7||lq , |T9|) 

ih 2 



(x,y\Q a (E) 
P*{E)\x,y) 



rn 



-K^ 2 (E)(G a (0, 0, E) - G fre e(0, 0, E))K l l 2 {E), 



n 



\/—im 

h 



9{ay){x\G a {y^E)K 1 ' 2 {E) 



e(ay)K 1 / 2 (E)G (T (Q,y,E)\x}, 



yj—im 

where the fourth propagator Go(y, y', E) is already defined in terms of the scattering 
resolvents ( |3~4| ) and 6(y) is the well known Heaviside step function. Then define a 
linear operator called half- derivative with the prescription 

d]j 2 e ay = a l/2 e a \ Rea 1/2 >0, (36) 

which is a sensibly defined positive square root of the differential operator d y . This is 
nonlocal operator which can be explicitly expressed for functions f(y) which increase 
slower that the square root as y goes towards plus or minus infinity, 



d l y /2 f(y) 



y , , f(y)-f(y') 

- {y-y'f 12 ' 



respectively. Thus for an operator valued exponential function, d^ 2 e lKy = \fiK 1 i 2 e lKv . 
One may use the Schrodinger equation with proper boundary conditions (which were 
used to derive eqs. (|l7|Jl9D ) to see that the scattering Green functions with one coor- 
dinate in the waveguide may be written as exponential functions in that coordinate 



G„(y, y >,E) = e-^yfW,E), 

<W, v, E) = g{y', E)e~^y, ay ~ U ' ay > U ' 



(37) 



for some operator valued functions / and g. Using the forms ([16|j37|) and the defi- 
nition of half- derivative one may rewrite the propagators in a more useful form 

-2 



{x\T a (E)\x'} 
(x,y\Q a (E)\x'} 

{x'\P a (E)\x,y) 



am 



d^d^ixMG^E) - G iTee (E)\x',y')\^, (38) 

(39) 



-d y i {x,y\G a (E)\x ,y }\ y >=o, 



y/orn 
M(ay) al 



df(x',y'\G a (E)\x,y)\ 



y'=0- 



(40) 



From these formulae one can easily derive semiclassical approximations by using the 
leading order semiclassical approximation for the energy dependent Green function 
(see e.g. [g, 0) 

9-7T 

£ Bj(q, q', E)e lS ^'' E ^ h ' w ^ 2 , (41) 



(q\G a (E)\q f ) 
B 3 (q,q\E) 
Dj(q,q',E) 



(2tt^)(/+ 1 )/ 2 

1 d E d ql S j dfeSj 



1/2 



m 



Ej(q, q',E) 



p yj p. 



//./ 



1/2 ' 



\detd x d x ,Sj(q,4,E)\ 



1/2 
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where the sum is taken over (usually finitely many) classical scattering trajectories 
labeled by j with classical actions Sj(q,q',E) = Jjdq ■ p, and Morse indices Vj 
which count the number of conjugated points along the orbit j. p y j = d y Sj and 
p' y j = —dy'Sj are the perpendicular (w.r.t. CSOS) projections of the final and 
initial momenta. Thus using the definition fl36|) in the leading semiclassical order 
the half derivatives only cancel the square roots of y-momenta if one expresses root 
of (/ + 1) x (/ + 1) determinant Bj in terms of root of (/ — 1) x (/ — 1) determinant 

Da 



{x\f a (E)\x'} = {2mh y f _ 1)/2 E D^x, 0, x', 0, E ySA^,o,E)/n^/2^ (42) 



[*,V\4*WW = ( 2^)// 2 g V* ^ ^^y^'\ (43) 



iy/2wme(ay) ^ Dj(x,y, x', 0, E) j S] ( x ,y, x ',o,E)/h-iu^/2 

Pyj 

WP„mx,y) S E ^^e«i^W, (44) 



(2mh)f/ 2 



In ( f4*2"|) the sum X)j is restricted only to classical orbits which strictly leave CSOS 
and lie entirely on the a-side of CS. The trivial classical scattering orbit whose 
y-coordinate is constantly zero is the only classical orbit of the semiclassical free 
Green function Gf Tee (E) and is thus canceled in the expression ( |38] ) for the propa- 
gator T a (E). The sums in (|43|j44] ) contain only the classical orbits which lie entirely 



on the a— side of the CS with one end point on the CSOS and the other end point 
lying in the a— side of CS. If there was a classical orbit whose part would lie in 
the waveguide then its y-coordinate should have an extremum there which, how- 
ever, is impossible since the classical motion in the waveguide is free in y-direction. 
The semiclassical expression for the CS-CS propagator without crossing the CSOS 
(x, y\ Go (E) | x', y') thus according to definition ([34]) looks the same as RHS of fl4lD 
where the sum now includes only the classical orbits which do not cross CSOS. The 
leading order asymptotic results ([|^,[|I||4]) agree with the semiclassical theory of 
Bogomolny [[j]]. 

The higher order semiclassical corrections to CSOS/CS-CSOS/CS propagators can 
be obtained in a systematic way by (i) inserting corrected higher order semiclassical 
expression for the scattering Green function (|4l| ) in the formulae ([3^, |39|, f|0D and 



(ii) evaluating the half-derivatives in terms of a power series in h. I will show now 
briefly how both steps can be systematically performed. 

(i) The higher order corrections to semiclassical energy- dependent Green function 
(q\G a (E)\q') can be obtained by multiplying each term of the formula ([|l]) by a 

5 A variant of this approach for time-dependent quantum propagator (without consideration of 
conjugated points - short time limit) has been developed by Roncadelli |l3|], whereas Gaspard and 
Alonso 0] used another (path integral) approach to derive ^.-expansion of Gutzwiller trace formula. 
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correction factor 

M G °WW = (2nih)U+W ^ B ^' h - 1 ^ 12 X^/fo, E). 

The corrections fl can be calculated by inserting the whole expression to the 
Schrodinger equation. Comparing the terms with equal power of h and integrating 
along the orbit yields the explicit recursion formulae for the semiclassical corrections 

fi( qj (t),4,E) = 1, 

fi( qj (t),q',E) = p| o t ^(-)^^(*')^^)A j/ ^ 1 (q j (t') ! g', J E), 

KjfM,E) = B- 1 (q,q',E)d q [B j (q,q',E)f(q,q',E)}, 

where qj{t) denotes the classical orbit j with end points q' and q. One must use the 
sign factors (— ) nUj and the Cauchy principal value of the integral in order to avoid 
infinite contributions each time one passes a singularity - conjugated point. 

(ii) Half derivative of a term like e tS ^ h f, where we shall take S = Sj, f = Bjfl, 
may be represented as a power series 

d^e^f = -±=e iS ' n £ h n C n f 
v n, n=0 

where C n are some linear operators independent of H. Taking it twice, d^ 2 d^ 2 = d y , 
and comparing the terms with the same power of H one obtains the set of equations 
which determine the operators C n 

Cq = idyS, 
CqC\ + C\Cq = d y , 

n 

C m C n - m = 0, n > 2. 

m=0 

It is easy to see that C n is n— th order differential operator. For example, we give 
explicit expressions for the first two 

C, = (45) 

and the next-to-leading order semiclassical expression for the quantum CSOS-CSOS 
propagator 

{x\f a {E)\x>} = f2 ^ 1)/2 E¥ SiW (46) 

3^,^- ifi^Dj /V7, ; /; ; V 

where all functions on RHS have arguments (as, 0, x', 0, E 1 ). 
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2.6 Symmetry of the CSOS-CSOS propagator 

At a given value of energy E one can split the SOS-Hilbert space on two orthogonal 
components, 

C= C Q (E)®C C (E), 

the (usually finite dimensional) subspace of open channels and the (usually infinitely 
dimensional) subspace of closed channels 

C Q (E) = {\i/>}e£, {m 2 (E)m>0}, 

c c {e) = {\ij}ec,{m 2 (E)m<o} 

spanned by {\n},E' n < E}, and by {\n},E' n > E}, respectively 

A very useful information about the CSOS-CSOS propagator which can be writ- 
ten in a block form 



± oo- 1 oo 


— T t f 

- 1 oo A oo 


if T 1 " 


— f 

- L oc ? 


7'T t f 


= f 

± CO ) 


if T t 


— ?T t f 



T co (E) T CC (E)J 

can be obtained by comparing the two expressions for scattering wavefunctions, 
conjugated eq. (§) and eq. (§). Comparing the values and normal derivatives on 
CSOS one obtains two equations 

VTiK Tl/2 (l±fa)\ip} = V±iK Tl/2 \l±flW}- 
Noting that K ±l / 2 K^ 1 / 2 ^ = f and performing some algebra yields 

1, (47) 

(48) 
(49) 

-fce-fl (50) 

This is the so called generalized unitarity [|TJ], [EJ of a CSOS-CSOS propagator. Note 
that the open-open part f°°(E) is indeed a unitary operator fl47j). 

Now we give the representation of the CSOS-CSOS propagator in terms of what 
is called reactance matrix in scattering theory ||. This gives also a practical recipe to 
determine the CSOS-CSOS propagators f a (E). Let ty an (x, y, E) denote the unique 
wavefunctions which satisfy Schrodinger equation (^) with boundary conditions 

ty an (x,0,E) = {x\n}, ^an(x, aoo, E) = 0, 

where the second condition should be taken on the boundary of CS if the latter is not 
infinite. Then one can define the reactance operators R a (E) with matrix elements 

{l\R„(E)\n} = — — ^— — - ( dx% l (x,y,E)d y ^ an (x,y,E)\ y=0 (51) 
V ( E ) kn ( E )s 

and show (using (^)) that the CSOS-CSOS propagators can be written as 

t a {E) = (1 + tR„(E))(l - iRa(E))- 1 . (52) 
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In case of time reversal symmetry the wave- functions ty an (x, y, E) are real and one 
can use Green's theorem to show that then the reactance matrix is symmetric 

{l\R a (E)\n} = {n\R a {E)\l} if <& an (x,y,E)=% n (x,y,E). 

Using representation (|52l) this means that then also the CSOS-CSOS propagator is 
symmetric 

{l\f a (E)\n} = {n\f a (E)\l}, or {x\f a (E)\x'} = {x'\f a (E)\x}. (53) 



2.7 Practical applications and semi-separable systems 

Let us truncate the basis of C to include all N Q open channels of C Q and the first iV c 
closed channels. The truncated (N = N + iV c )— dim matrices with matrix elements 
{l\T a (E)\n} and {l\R a (E)\n} will be denoted by T a (E) and R a (E), respectively. In 
practice one should increase N c until numerical results converge, which is typically 



the case [14, O, [TT| for very small values of N c already, so in the semiclassical limit 



iV « N . The practical SOS-quantization condition then reads 

det(l-T i (E)T T (E)) = 0. (54) 

Using ([52]) this condition can be formulated in terms of reactance matrices 

det(R T (£) + Ri(E)) = 0. (55) 

In case of systems having a time-reversal symmetry (that is if H'(p x ,x,y) = 
H'{ —p x , x, y)) the reactance matrix ( pi]) is a complex symmetric matrix R a (E) = 
R^(E) which is written in terms of a purely real matrix ~R a (E) via 



where the diagonal elements are N and A^ c dimensional sub-matrices. The quantiza- 
tion condition (|55| ) can be thus expressed in terms of purely real symmetric matrices 

R a (E) = R T a (E) = Rl(E) 

det(R T ( J E) + Ri(E)) = 0. (56) 



The equation (|56"D is much more efficient for numerical calculation of energy spectra 



(by seeking its zeros) than original quantization condition fl54|) since the former 
involves real arithmetic ||11|| . 

There is a generic (in a sense of dynamics) class of the so called semi- separable 
systems for which reactance matrices can be straightforwardly calculated and hence 
quantization condition (|56|) can be easily implemented. Semi-separable system 
should be separable (in (x, y) coordinates) on both sides of CSOS (for y > and 
y < 0) but it can be discontinuous on CSOS (H'(— 0) ^ H'(+0)). Thus one has two 
complete sets of normalized SOS eigenmodes, first \n} + are eigenstates of H'(+0), 
and the second |n}_ are eigenstates of H'(— 0). Since the system is separable on 
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both sides one can explicitly calculate the wavefunctions *S? an (q, E) by separation of 
coordinates 

^ 1n (x,y,E) = (f) 1n (y,E){x\n} + , y > 0, 
^ln{x,y,E) = <p ln (y,E){x\n}-, y < 0, 

where y-dependent parts should be normalized to give 4> an (0) = 1. We have free- 
dom to cut CS slightly above discontinuity and choose priveliged set \n} + with 
wavenumbers k n (E). Then we apply the definition (ETT) to calculate real reactance 
matrices 

R ul {E) = \k n {E)\- 1 d y ^ n {0,E)5 nl , 

R lnl (E) = —\k n (E)ki(E)\~ 1 ^ 2 ^2 +{ n \J}- 9y4>u(0, E) 

j 

The upper is diagonal while the lower includes transformations by means of real 
orthogonal matrix with elements _{Z|n} + = which are typically easily 

calculated knowing the two sets of SOS-eigenmodes. The author has applied this 
method for numerical calculation of energy levels in the so called semi- separable 2- 
dim oscillator |ll| which is a generic nonlinear autonomous dynamical system with 



two freedoms. The method turned out to be capable of yielding consecutive energy 
levels with sequential number of the order of 10 6 (million !) at the cost of few 
minutes of Convex C3860 CPU time per level. 



3 Abstract formulation of the method 

In this section we devise a general and abstract mathematical framework within 
which one can prove all versions of simply connected SOS quantization conditions 
and SOS decompositions of the resolvents of the corresponding Hamilton operators. 

Let M. be an arbitrary normed vector space, which will be referred to as reduced 
space. The vectors from reduced space Ai - r-vectors will be written in bold and 
linear operators over reduced space - r-operators will wear mathematical accent ". 
Then we define r-operator valued function h(y, E), where y is a real variable and E 
is a complex spectral parameter (e.g. energy), in order to study the following general 
homogeneous vector differential equation 

(d y -L(y,E))f(y) = (57) 

on the entire real axis y G TZ (or on some finite interval which contains zero). 
Normalized r-vector valued functions f(y), J dy\\f(y)\\ 2 < oo, constitute a normed 
vector space, denoted by TL. The values of spectral parameter E for which the equa- 
tion (|57D has nontrivial solutions in 7i are called generalized eigenvalues whereas the 
corresponding r-vector valued functions are called generalized eigenfunctions of the 
equation flS?]). It will be shown in the next subsection that this problem is equivalent 
to time-independent Schrodinger equation for some special choice of L(y,E). The 
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equation ([57]) can be generally solved by means of a Green function G(y, y', E) which 
is a unique r-operator valued function which solves the inhomogeneous equation 

(d y - L(y, E)) G(y, y') = 5(y - y')J (58) 

with boundary conditions 

lim G(y,y',E) = 0. (59) 

y — >±oo 

where J can be an arbitrary nonsingular r-operator. 

Equation ([58]) may be written in the form L(E)G(E) = 1 where L(E) and 
G(E) are the operators over 7i with kernels being r-operator valued distributions 
J (5'(y—y')—~L(y,E)5(y—y')) and G(y,y', E), respectively. If left and right inverse 
of L(E) exist and coincide then G(E)L(E) = 1, so the Green function satisfies also 
the "conjugated" equation 

G(y, y') (d yl - L V, E)) = -5(y - y')J. (60) 

where 

i\y,E) = -r 1 i(y,E)3. 

We shall construct the Green function G(y,y', E) by means of the Green func- 
tions G a (y, y', E) of two generalized scattering problems which are defined by cut- 
ting the y— axis at y = and substituting the upper (y > 0, o =|= +)/lower 
(y < 0,cr =[= — ) part of function L(y, E) by a constant h(E) = L(0, E). Therefore 
these scattering Green functions satisfy 



d y -Hy,E))G.{y,y') = S( y - y ')J, if ay > o, (61) 

(d y -L(0,E))G a (y,y') = S(y-y')J, if ay < 0, (62) 



with boundary conditions 



lim G a (y,y',E)=0. (63) 

y— >ctoo 



The scattering Green function G a (y, y', E) can be written explicitly on (—a)— side 
{ay < 0, ay' < 0) in terms of the abstract scattering operator T cr (i?) 



G a {y,y',E) =tex V (t(E)y) (| - i[y - y'}J - T a (E)) exp(L (E)y'), (64) 

where [y] =t= +]y > 0, [y] =|= — ;y < denotes the side or sign. The equation 
( |64"D can be considered also as a unique definition of the abstract scattering operator, 
or 

T ff (E)=iG a (±0,^0,E) + ^iJ. (65) 
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Theorem 2b: The Green function of the equation $58 ) G(y, y',E) can be decom- 
posed in terms of four r-operator valued functions which can be defined by means of 
the scattering Green functions 



G (y,y',E) 

P(y,E) 

T(E) 



G[y](y,y',E), 



(66) 

Vi[y]G [y] (y,0,E), (67) 
Vi[y]G [y] (0,y,E), (68) 
%(E) + t^E) = *G T (±0, =f0, E) + zG^TO, ±0, E) + 1. (69) 

Namely, the decomposition formula reads 



G(y, y', E) = G (y, y', E) + Q(y, E)(l- T(E) P(y', E) 



-l 



(70) 



Proof: One must show that RHS of Q7DD is the solution of equations ( j5q , p9| ). The 



first term of RHS is indeed a solution of a inhomogeneous equation fl58|) with bound- 
ary conditions Q59D on both sides but it is generally discontinuous at y = 0. The 
second term of RHS, or its first factor Q(y,E), is a solution of the homogeneous 
equation (|57D on both sides but it is again discontinuous at y — 0. The sum of the 
two terms is therefore also the solution of the inhomogeneous equation (0). One is 
left to prove that the sum is continuous at y = and therefore a unique solution of 
equation (0). For arbitrary function of y we define the difference operator 

\f(y) = f(y + 0)-f(y-0). (71) 

Then the straightforward calculation yields 

A y G(0,y',E) = A y G (0,y',E) + (A y Q(0,E))(l-T(E))- 1 P(y',E) = 
= A y G (0,y',E) - [y']G lyl] (0,y',E) = 0, 

since 

Aj,Q(0, E) = V~i(l - T(E)). (72) 



Theorem lb: For any generalized eigenvalue E of the equation ( \5l\ ) the opera- 
tor T(E) has a fixed point. 

Proof: Any generalized eigenfunction of the equation (|57D corresponding to gener- 
alized eigenvalue E can be written in a form 

f(y) = Q(y,E)a (73) 

for some nonzero r-vector a E Ai. One can write explicitly, a = Q 1 (0,-E')/(0) if 
Q(0, E) is invertible, or more generally, a = [d r y Q(y, E)]- l d r y f(y)\ y=0 if <9£Q(y, E)\ y=0 
is singular for all p = 0, 1 . . . r — 1. Eq. ([73) follows from the definition of r-operator 
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valued function Q(y,E) in terms of scattering Green functions on the nontrivial 
side. Since the function f(y) is continuous at y — 0, A y f(0) = 0, one sees, using 
equation (|7^), that a is a fixed point of the operator T(E), 

T(E)a = a. 

Note that general decomposition formula ( |7D| ) is invariant with respect to simi- 
larity trasnformations 

Q(y,E) - Q(y,E)S, 

T(E) -> S X T(^)S, (74) 

and transformations 

Q(y,E) - Q(y,E)Zr\ 
T(E)-1 -> Z(T(£)-1)Z, (75) 

P(y,£) - ZT'p^E) 
where S, Z are any bijective r-operators. 



3.1 Trivial application 



For example, let us first cast our ordinary Schrodinger problem (§,@) into the abstract 
form. Reduced space should be here M. = C © £, since the Schrodinger equation is 
of the second order. One should take 



-l 





K 2 (y,E) 



-1 




(76) 



in (|3J5|) where K 2 {y,E) = {2m/h 2 ){E - H'(y)), K{E) = K(0,E). Then, refer- 
ring to the two components C of M. with indices 1 and 2, \i/j(y)} = fx(y) is a 
solution of the Schrodinger equation (d 2 + K 2 (y, E))\ip(y)} = and G(y,y',E) = 
— (2m/h 2 )Gii(y, y\ E)K~ 1 (E) is its normal Green function (13) in hybrid repre- 
sentation. Comparing scattering Ansatze (|T3) and ( Ell ) and using the similarity 
transformation (ff^) 

■R-^iE) k~^ 2 (E) 
l K 1 / 2 (E) -iK l / 2 (E) 

one obtains 

One can also check that the first "row" of (ih/ \ / 2m)Q(y, E)S and the first "column" 

of (^/ v / 2^)S _1 P(y, E)K(E) can be written as (Q^E), Q T (£)) and (P T (£), P;(£)) T , 
respectively. It is now easy to check that 1, 1 component of the general decomposi- 
tion formula ([70]) agrees with the more special case (|35|). 
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3.2 Nontrivial applications 



There is a straightforward nontrivial generalization of application (|76|) , namely, one 
can include nonrelativistic systems which interact with very general external (gauge) 
fields and have thus Hamiltonians in our canonical coordinates (x, y) of the form 

H = ^HH - Mv)) 2 + H'(y), (78) 

where the only restriction for the self-adjoint operators A(y) = A(—ihdx, x, y), H'(y) 
H'(—ihdx,x,y) is that they should not depend upon p y so they can be restricted 
to act over small SOS-Hilbert space C. Again we define the restricted space as 
Ai = C © C and Schrodinger equation (|j) with (|78| ) can be written as a first order 
system (fTj) where (|76|) should be replaced by 



L ^^-{K^y,E) iA(y)/h)> J -{l 

where all statements from the previous example remain valid except that now the 
CSOS-CSOS propagator cannot be separated to upper and lower part like (|77|) and 
all blocks of T(E) are generally nonzero. 

As for another interesting application one can decompose the Green function of 
a relativistic Dirac |— spin fermion bound in external electromagnetic field A^(x,y) 
and search for its stationary states. One may choose e.g. x = (x 1 ,^ 2 ),?/ = x 3 and 
L(y, E) = -f 3 (-f°(-iE + ieA )+-f 1 (d 1 -ieA 1 )+^ 2 (d 2 -ieA 2 ) + im) + ieA 3 J J = -7 s , 
so that equations (|57| , |58|) reduce to Dirac equation where % — c — 1. The reduced 
Hilbert space is now the space of Dirac spinor valued functions over 2-dim plane 
(x 1 , x 2 ). 



4 Multiple sections 

In this section I consider the case of multiply connected CSOS. Let CSOS, which 
is now a smooth multi-sheeted (/ — l)-dim manifold S, divide the /—dim CS C on 
countably many disconnected parts whose closures are denoted by C a , a G J , 

C= \JC a 

where J is some finite or countable index set. Two points are in the same compart- 
ment C a if they can be connected by a continuous curve which does not cross CSOS 
S. The compartments C a and Cp are said to be neighbouring (denoted by a\(3) if 
their intersection S a p is nonempty (/ — l)-dim manifold 

a\P & C a nCp = S aP + 0. 

The union of all such intersections is the whole CSOS 
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H is self-adjoint operator over the Hilbert space H. = L 2 (C). Let O a p, a\/3 be such 
open sets which cover the connected parts of CSOS, S a p C O a p. Hamiltonian 
operator TC is admissible if there exist coordinates (x,y) a p for each of the sets O a p 
such that 

H\mo afj ) = -\fi 2 d y M^{y)d y + H' ap {y), (79) 
M a p{y) = M Q p(-ihd x ,x,y), 
H' a p(y) = H' a/3 (-ihd x ,x,y). 

We choose the sign of coordinate y of O a p so that (x,y) a p G Cp if y > 0. Here we 
have allowed for very general "masses" M a p(y), which should of course be positive 
operators and hence invertible, which is another generalization of this section. Then 
I introduce small SOS-Hilbert spaces C a p = L 2 (S a p), a\(3. The operators restricted 
to C a p will be again denoted by accent Now cut off the CS around C a and attach 
y— flat the so called a/3-waveguides on the other sides of all connected parts S a p, a\(3 
of the boundary dC a (see figure 2). Thus one defines the scattering Hamiltonians 
which in local coordinates read 

^I^)-(_ (r72H m-^oh + ^(o); W ' l ) 

The fundamental solution of the time-independent Schrodinger equation for the 
scattering problem (^) in the a/3-waveguide is given by 

,Jx|M^ 2 (0)^ 1/2 (E)e ± ^(^ 

where the wavenumber operator K a p(E) is the positive square root of the self adjoint 
operator 

KlpiE) = |M^ 2 (0)(i? - ^(0))M^ 2 (0). (81) 

Vectors from the dual space £' a/3 are written with reversed indices, e.g. p a {x\ E 
C a p- Thus the general scattering wavefunction of the Hamiltonian H a in the a/3- 
waveguide (y > 0) reads 



-jf p a {x\M l Jp 2 K-T' 
+ e ik ^ E ^fp ai {Em a ± (82) 



7 



K(x,y,Er) = ^ 1 -[p a {r\e- l ^ {E)y + 



+ E 7a {^lT 7 ^(E)e^^]^ /2 (E)M^ 2 |^} Q , ; (83) 

7 

(where M^p 2 = M^ 2 (0)) and is uniquely determined by the incoming waves para- 
metrized by the SOS-states |^} a7 or l^*}^ coming from the «7-waveguide for all 
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neighbouring compartments C 7 . We have introduced the scattering operators T]g Q7 
which will be called generalized CSOS-CSOS propagators. Tp ai is the scattering 
operator from £ Q7 to £p a and describes the propagation from C 7 to Cp via C a . Then 
we define the two types of linear operators: Q' ai from small SOS-Hilbert spaces C ai 
to Hilbert space TC, and P^ from Hilbert space Ji to small SOS-Hilbert spaces C ia 
by the following prescriptions 

The resolvent of the scattering Hamiltonian with outgoing boundary conditions 

G a (E) = {E-H a + i^)~ 1 



can also be written explicitly (in analogy with equation fllq)) inside the waveguides 
(y>0,y'>0) 

((x,y) a p\G a (E)\(x',y') ai ) = i fr{x\ti$&$*(Ef)[6^>*W\y-*\ + 

in 

+ e^^% a7 (£)e^^ (84) 

In analogy with the simply connected case we define also: (i) The operators Q aj (E) 
from C ai to H and the operators P ya (E) from Ti. to C ia with the kernels 

(<7l&* 7 (£)IV}a 7 = JtelQ^WW; <?ec a , (g5) 

[ U, q £ u a , 

^1^(^19) = q q f- (86) 

which are called generalized CSOS-CS and CS-CSOS propagators respectively, and 
(ii) Gq{E) a linear operator over 7i with the kernel 

cmw) - { 3 2*t e c " ^ 

which is called generalized CS-CS propagator (without crossing the CSOS in be- 
tween) . 

Let us compactify our notation by introducing the following symbols. First we 
define the large SOS-Hilbert space M. 

M=($£ af3 
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with a complete system of orthogonal projectors Tl a p- (Note that each pair (a, ft) 
is always included twice, once as a\(3 and once as f3\a.) For each SOS-state \ip} we 
write symbolically 

1^} = I^W, |^}a/3 = tl a/3 \lp}. 

One then defines the large operators f(E), Q(E) and P(E) by 

T{E) = Y Tp ai (E)Tl ai , 

/3|a| 7 

q{e) = £Q Q7 (p)n a7 , 
p(p) = ]Tn 7a p 7Q (P), 

or equivalently 

Up a f(E) = Y.f Pai {E)ti ai) (88) 
Q{E)ti ai = Q a7 (P)fl a7 , (89) 

n 7a p(P) = n 7a p 7a (P). (90) 

Now, geometrically most general form of the main result of this paper can be stated 
and proved very elegantly. 

Theorem 2c: The resolvent of the Hamiltonian G(E) — (E — H)^ 1 can be decom- 
posed in terms of four elementary propagators, namely G : H — > Ti, Q : M. — > H, 
Q : H -> M, and f : M -> A4, as follows 

G(P) = G (P) + Q(P) (l - T(P)) ~* P(P). (91) 



Proof: Put the decomposition formula in a sandwitch between (q\ and One 
should prove that RHS also solves the inhomogeneous Schrodinger equation as the 
LHS does 

(P - H(-ihd q , q))(q\G(E)\q') = S(q - q'). (92) 

This is indeed true in every compartment C a separately, where (q\Go(E)\q') is a 
particular solution and (q\Q(E)(l — T(E))~ 1 P(E)\q') is a solution of the homoge- 
neous equation by construction of operator Q(E). What is left to prove is that RHS 
is continuously differentiable on borders between compartments, that is on CSOS 
S. Take arbitrary neighbouring compartments C a and Cp and choose coordinates 
(x,y) a p of an open set O a p which includes S a p. We shall need the following values 
and normal derivatives of the CSOS-CS and CS-CSOS propagators on the S a p which 
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can be obtained directly from (|^,|3|,[84], [84 



((SB, 0)up\Qary(E) 

dy{(x,0) a p\Q aJ (E)\ y=0 
P ia (E)\(x,0) aP ) 
d y P ia (E)\(x,y) a p)\ y =o 



/3a 



{x\M l J p 2 K-l ,2 {E) (f Pa7 (E)+5, 



Pi 



(93) 
(94) 

X^ + d^K^WM^lx}^, (95) 



\ (f jaP (E) - Vy) Kl^EW^lxU. (96) 



First we shall prove that RHS of (|9~1~D is continuous on all S a/ 3, a\(3. Using the 
difference operator (|7T| ) we can write (omitting the argument (E) for the sake of 
brevity) 

A y ((x,0U\G\q') - A y ((x,0U\G \q') = 

(7|a -y 1/3 \ 

£ Q ai ti aj - E 4^ (i - fy'Piq') = 



1 fnlM 1 / 2 ^- 1 / 2 

f3a 



{x\M^K~p 



2J (^/3<*7 + ^/3 7 )n«7 — E (^a/37 + ^a7)fl/3 7 
7GJ 76^ 



(l-T)-^!^) 



^{x\M^ 2 k-; /2 \P af3 -Pp a \\q'). 



We have applied equations (|93| , |88| , |90f ) . And analogously, by applying equations 
( p4] , p8| , pOD we get (note also that (x,y) a/3 = (x, —y)p a ) 



Ayd y ((x,y) al3 \G\q')\ y= o - A y d y {(x,y) a/3 \G \q')\y =0 = 

y\a ^ ^ j\/3 

2J (^3«7 ~~ ^/3 7 )ria7 + 2J (^Qi/37 — ^7)^/37 
7G^ 76^7 



(l-TT^FV) 



% r , ,- f l/2 t>1/2 



ft 



-fa/3 + Pf3a 



\q')- 



In order to see that A y {(x, 0) ai g|G|q') = and A y d y {(x, y) a p\G\q') \ y= o = one has 
to prove 



<(x,0) a/3 |G Q (£)|< Z , > 
^<(a?,2/)oi9|G ! a(S)|g'>| w =o 



ft 
ft 



1 , a {x\M% 2 k-l'\E)P, a {E)\q'), (97) 



^{xiM^k^^P^mq'), 



(98) 



considering the definition of G (E) in terms of G a (E) fl8~T|). But this is easy Both, 
LHSs and RHSs of (|97lj98| ) satisfy the conjugated Schrodinger equation as functions 
of q'. The initial data — the values and the normal derivatives of LHSs and RHSs on 
any initial surface (x, y) ai , also match as can be seen by applying q' = (x,y) ary 
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and flU) to LHSs and (|95|j96|) to RHSs. Formulas (|3J9|)/(|T]) then follow from the 
uniqueness of the initial value homogeneus (f|) /non- homogeneous ( |92"|) Schrodinger 
problem, q.e.d. 

One can formally expand decomposition formula ( |91~1) in a geometric series or 
sum over paths 



n>2 



G{E) = G Q (E) + £ Q (S)T Q (E)P a2ai (E) 



a.x\a2—\un 

(99) 

where each term contains probability amplitudes to propagate from compartment 
C ai to C a2 ... to C ajv . If one chooses many disconnected parts of CSOS S ak0lk+1 which 
are geometrically close then the propagators T a/ g 7 (i?) would become simple and they 
could be asymptotically explicitly calculated, so the formula (|99|) would turn into a 
kind of path integral formula for energy dependent quantum propagator. 



5 Discussion and conclusions 

This paper presents theoretical construction of SOS reduction of quantum dynamics 
in analogy with the SOS reduction of classical dynamics. However, there is an im- 
portant difference: In classical dynamics, one should carefully choose SOS such that 
almost every trajectory crosses it, while in quantum dynamics this is not essential. 
All theorems work even if CSOS lies in classically forbidden region although the 
practical usefulness of the method is expected to be worse then, because of the ex- 
ponential localization and sensitive dependence on boundary conditions. Moreover, 
formalism of section 4 can be easily adapted (by taking two different CSOSs as a 
single multiply connected CSOS) to show explicitly that the spectra as determined 
by our method do not depend on the choice of the CSOS, since the correspond- 
ing quantum Poincare mappings are related to each other by a kind of similarity 
transformation. 

The Green function - energy dependent quantum propagator has been decom- 
posed in terms of propagators which propagate from CS/CSOS to CS/CSOS. This 
decomposition formula has been generalized in two ways: (i) for Green functions of 
arbitrary linear differential systems and (ii) for SOS which consists of more than one 
disconnected part. The combination of these two generalizations is straigtforward 
so it is not given explicitly in this paper. While this general decomposition formula 
((|9T|) or even (|70|)) has so far merely theoretical value, it has a very practical con- 
sequence, namely, SOS-quantization condition: The resolvent of the Hamiltonian 
{E — H)^ 1 can have a pole, i.e. eigenenergy E Q , only when the operator 1 — T(E ) 
is singular, i.e. when general quantum Poincare mapping T(E ) has a fixed point 
l^} G Ai, T(E )\ip} = For the more special and common case of section 2 we 

have M. = C ® T = "j)' ) ' an< ^ ^ = ( j || ) ' w ^ ere I T } is a fixed point of 

quantum Poincare mapping Tj.ff and at the same time | |} is a fixed point of a sim- 
ilar mapping IfT;. This quantization condition can be very efficiently numerically 
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implemented fL4 , [12|, [TT|]. Since the exact quantum Poincare mapping is usually 
difficult to calculate explicitly we describe its semiclassical ft— expansion and give 
explicitly the leading (Bogomolny's and next-to- leading order terms. 



Acknowledgments 

I am grateful to Professor Marko Robnik for fruitful discussions. The financial 
support by the Ministry of Science and Technology of the Republic of Slovenia is 
gratefully acknowledged. 



26 



References 

[I] Bogomolny E B 1992 Nonlinearity 5 805; 1990 Comments At. Mol. Phys. 
25 67 

[2] Doron E and Smilansky U 1992 Nonlinearity 5 1055 

[3] Feynman R P and Hibbs A R "Quantum Mechanics and Path Integrals" 
(McGraw-Hill, New York, 1965) 

[4] Gaspard P and Alonso D 1993 Phys. Rev. A 47 R3468 

[5] Lichtenberg A J and Lieberman M A "Regular and Stochastic Motion" 
( Springer- Verlag, New York, 1983) 

[6] Littlejohn R G 1990 J. Math. Phys. 31 2952 

[7] Ozorio de Almeida A M 1994 J. Phys. A: Math. Gen. 27 2891 

[8] NEWTON R G "Scattering Theory of Waves and Particles" (Springer- Verlag, 
New York, 1982) 

[9] Prosen T 1994 "Quantum Surface of Section Method: Decomposition of the 
Resolvent (E - H)^ 1 " Preprint CAMTP/94-3 submitted to J. Math. Phys. 

[10] Prosen T 1994 J. Phys. A: Math. Gen. 27 L709 

[II] Prosen T 1994 will be submitted to J. Phys. A: Math. Gen. 

[12] Rouvinez C and Smilansky U 1994 "A Scattering Approach to the Quan- 
tization of Hamiltonians in 2 Dimensions - Application to the Wedge Billiard" 
Preprint IPT Lausanne 

[13] Roncadelli M 1994 Phys. Rev. Lett. 72 1145 

[14] Schanz H and Smilansky U 1993 "Quantization of Sinai's Billiard - a 
Scattering Approach" Preprint Weizmann Institute 

[15] Weidenmuller H A 1964 Ann. Phys. (NY.) 28 60 



27 



Figures 



Figure 1 The geometry of 2-dim CS of a typical bound system (a) with simply con- 
nected CSOS. Isopotential contours are shown. The product of classical or quantal 
scattering mappings of the two scattering systems shown in (b,c) is equal to the 
classical or quantal Poincare mapping of a bound system (a). 

Figure 2 The geometry of 2-dim CS of a bound system with multiple sections (a). 
One of the related scattering systems is shown schematically in (b). 
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